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ABSTRACT 

We present a new catalogue of 5,106 infrared bubbles created through visual classi- 
fication via the online citizen science website 'The Milky Way Project'. Bubbles in the 
new catalogue have been independently measured by at least 5 individuals, producing 
consensus parameters for their position, radius, thickness, eccentricity and position 
angle. Citizen scientists - volunteers recruited online and taking part in this research 
- have independently rediscovered the locations of at least 86% of three widely-used 
catalogues of bubbles and H II regions whilst finding an order of magnitude more 
objects. 29% of the Milky Way Project catalogue bubbles lie on the rim of a larger 
bubble, or have smaller bubbles located within them, opening up the possibility of 
better statistical studies of triggered star formation. Also outlined is the creation of 
a 'heat map' of star-formation activity in the Galactic plane. This online resource 
provides a crowd-sourced map of bubbles and arcs in the Milky Way, and will enable 
better statistical analysis of Galactic star-formation sites. 

Key words: H II regions - infrared: ISM - ISM: dust - stars: formation 



1 INTRODUCTION 

H II regions ionised by young, O and B-type stars pro- 
vide the most readily observable tracers of star formation 
in the Milky Way and other galaxies. Ionised gas produces 
strong emission in optical and infrared (IR) recombination 
lines, forbidden lines and thermal (free-free) radio contin- 
uum. Dust mixed with ionised gas and heated by the hard 
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radiation field makes H II regions bright sources of thermal 
IR emission. 

In the Milky Way, the spatial morphology of individ- 
ual H II regions can generally be resolved, revealing complex 
structures which are often shaped by the newly formed stars 
(Anderson et al. 2011). One particularly common and infor- 
mative morphology is the presence of a rounded 'bubble' of 
emission from gas excited presumably by a central source. 

Churchwell et al. (2006, 2007, hereafter CP06, CWP07) 
catalogued nearly 600 IR 'bubbles' in the inner 130° of the 
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Galactic plane by visually searching for ring-shaped struc- 
tures (complete or partially broken) in 3.6-8.0 /im images 
from the Galactic Legacy Infrared Survey Extraordinaire 
(GLIMPSE; Benjamin et al. 2003; Churchwell et al. 2009). 
In the GLIMPSE-I survey area, CP06 catalogued 322 bub- 
bles, of which 25% corresponded to previously catalogued 
radio H II regions and 75% were attributed to late B-type 
stars with insufficient ionising luminosity to produce radio- 
bright H II regions. Bania et al. (2010) later selected 24 /im 
diffuse emission sources from Spitzer/MIPSGAL survey im- 
ages (Carey et al. 2009) with spatially coincident 21 cm con- 
tinuum emission for a radio recombination line survey with 
the Green Bank Telescope, and discovered several hundred 
new Galactic H II regions in the region —16° < / < 67° and 
|6| < 1°, including those associated with 65 of the CP06 
bubbles. 

GLIMPSE images are particularly helpful in identify- 
ing bubbles as emission from polycyclic aromatic hydrocar- 
bons (PAHs) at 8 /im demarcates the bubble rims in Spitzer 
Space Telescope images, while a second peak, arc, or torus of 
24 fim emission from warm dust is frequently observed inside 
the bubble with a morphology that closely traces the radio 
continuum emission (CP06; Watson et al. 2008, 2010). This 
qualitative pattern of bright 8 /J,m emission shells surround- 
ing regions of bright 24 /im/radio emission is also charac- 
teristic of giant Galactic H II regions and extragalactic star- 
forming regions (e.g. Povich et al. 2007; Bendo et al. 2008; 
Relano & Kennicutt 2009; Flagey et al. 2011). 

Nevertheless, existing surveys are not sufficient to es- 
tablish the true relationship between bubbles and HII re- 
gions. Both CP06 and CWP07 stressed that their bubble 
catalogues were far from complete, particularly with regard 
to bubbles with large (>10') and small (<2') angular di- 
ameters, and the majority of previously catalogued, bright 
radio H II regions (Paladini et al. 2003) in the GLIMPSE 
survey area were not associated with catalogued bubbles. 
None of the large bubbles identified by Rahman & Murray 
(2010) were included in the CP06 and CWP07 catalogues, 
nor was the large bubble associated with the well-studied 
giant H II region M17 identified by Povich et al. (2009). 
These large bubbles were missing from the catalogues be- 
cause they are generally faint, broken, and confused with 
smaller, brighter H II regions. 

In this paper, we present a new catalogue of 5,106 IR 
bubbles identified via visual inspection of the GLIMPSE and 
MIPSGAL survey images as part of the citizen science Milky 
Way Project 1 (MWP). Our catalogue expands the CP06 
and CPW07 catalogues by nearly an order of magnitude 
and hence represents a far more complete sample of Galactic 
H II regions. Three key advances in the bubble selection pro- 
cess enabled the identification of bubbles that were missed 
by CP06 or CWP07: (1) We enlisted >35, 000 volunteers to 
examine the images rather than relying on only a handful of 
experts; (2) we incorporated MIPSGAL 24 /im data, which 
greatly facilitates the identification of bubbles compared to 
relying primarily on 8.0 fim data; and (3) we relaxed the cri- 
teria of what defines a 'bubble' to include more incomplete 
shells and arc-shaped structures. The remainder of this pa- 
per is organised as follows: In sections 2 and 3 we describe 
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Figure 1. Screonshot of the Milky Way Project user interface. 
Colour figure available online. 



the MWP and the construction of our catalogue. Some ini- 
tial results from the catalogue are presented in section 4 
and discussed in section 5. We summarise our conclusions 
and directions for future work in section 6. 



2 THE MILKY WAY PROJECT 

Galaxy Zoo (Lintott et al. 2008, 2011) and the larger suite 
of Zooniverse 2 projects (Smith et al. 2011), have successfully 
built a large community of volunteers 3 eager to participate 
in scientific activities. 

The Zooniverse has shown that enlisting 'citizen scien- 
tists' via the internet is a powerful way to analyze large 
amounts of data. Human brains excel at pattern recogni- 
tion tasks, and most people will reach a level of accuracy as 
high as any expert after a brief introduction. By enlisting 
citizen scientists, researchers can extend visual classification 
to large samples of images, having each image examined by 
a large number of independent classifiers. This allows re- 
searchers to tap into the 'wisdom of the crowd' effect where 
the consensus of a group of non-experts is often more accu- 
rate than the testimony of a single expert. 

Another advantage of enlisting human classifiers is their 
ability to recognize unusual objects which computer search 
algorithms may be unable to spot. This has been shown by 
the serendipitous discovery of Hanny's Voorwerp (Lintott 
et al. 2009) and the case of the Galaxy Zoo 'Green Peas' 
(Cardamone et al. 2009). Spitzer GLIMPSE data is ideally 
suited to classification by citizen scientists as the amount 
of data is large and the images contain complex, overlap- 
ping structures that are impossible to disentangle using au- 
tomated algorithms. The task of recognising bubbles may 
eventually be handled by advanced machine-learning algo- 
rithms (e.g. Beaumont et al. 2011) but in the meantime the 
community of Zooniverse users are keen to contribute to as- 
tronomy and science (Raddick et al. 2010). 

The MWP is the ninth online citizen science project 
created using the Zooniverse Application Programming In- 
terface (API) toolset. The Zooniverse API is the core soft- 
ware supporting the activities of all Zooniverse citizen sci- 
ence projects. Built originally for Galaxy Zoo 2, the software 

2 http://www.zooniverse.org 

3 Over 480,000 registered volunteers at time of writing. 
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is now being used by 11 different projects. The Zooni verse 
API is designed primarily as a tool for serving up a large 
collection of 'assets' (for example, images or video) to an in- 
terface, and collecting back user-generated interactions with 
these assets. 

The assets in the MWP are multiband, false-colour 
JPEG images, created by gridding the Spitzer GLIMPSE 
and MIPSGAL mosaics into smaller images at three dif- 
ferent zoom levels. The highest zoom level provides users 
with tiles of 0.3° x 0.15°, and at a resolution of 800x400 pix- 
els these tiles nearly reproduce the 1.2" pixel scale of the 
GLIMPSE survey images. Larger tile sizes of 0.75° x 0.375° 
and 1.5° x 0.75° were also generated. The tiles were plotted 
in an overlapping grid to allow all parts of the inner Galactic 
plane (| l\ < 65° , |6| < 1° ) to be viewed by the MWP users, at 
all zoom levels. To provide an optimal representation of the 
dynamic range within each tile, each of 3 single-band images 
was independently scaled to a square-root stretch function 
(with the faintest 5% of image pixels clipped to black and the 
brightest 0.2% clipped to white), assigned to a colour chan- 
nel (red=24 /im, green=8.0 /jm, blue=4.5 (im), and finally 
composited into a 3-colour image. The MIPSGAL 24 /itm 
mosaics frequently saturate in regions of bright nebulosity, 
and saturated 24 /j,m pixels were set to maximum red to 
preserve the visual appeal of the images and to avoid pre- 
senting MWP users with saturation artefacts. The resulting 
composite images allow visual identification of both bright 
and faint features within a given image tile. 

The MWP user interface (see Figure 1) was built using 
Flash, based upon the pre-existing Moon Zoo interface (Joy 
ct al. 2011). Volunteers are primarily encouraged to draw 
ellipses onto the image to mark the locations of bubbles. A 
short, online tutorial shows how to use the tool, and exam- 
ples of prominent bubbles are given 4 . As a secondary task, 
users can also mark rectangular areas of interest, which can 
be labelled as small bubbles, green knots, dark nebulae, star 
clusters, galaxies, fuzzy red objects or 'other'. Examples of 
these are also given in a tutorial on the website, and these 
are discussed further in Section 3.4. Users can add as many 
annotations as they wish before submitting the image, at 
which point they are given another image for annotation. 
Each image's annotations are stored in a database as a clas- 
sification, and users can see the images they have classified 
in a part of the site called 'My Galaxy'. Users can only clas- 
sify a given image once. 

When marking bubbles, users place a circular annulus 
that can be scaled in size and stretched into an elliptical 
annulus. As they first draw out an object, they are able to 
control the position and size of the bubble. Once the bubble 
has been drawn they can edit these initial parameters as well 
as the bubble's ellipticity, annular thickness and rotation. 
In this way users can attempt to match the bubbles they 
see in each image to give an accurate representation. Users 
are also able to mark regions of the annulus where there is 
no obvious emission - as in the case of broken or partial 
bubbles. These 'cut-outs' are created by erasing (and then 
re-filling if necessary) segments of the annulus. Of all the 
bubbles drawn, 75% have thicknesses other than the default 
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Figure 2. Distribution of users with number of classifications. 



or minimum values, 50% are non-circular (of which 56% have 
been rotated) and only 12% were drawn with cut-outs. 

When marking another area of interest on an image (e.g. 
star clusters, green knots, etc) users simply draw rectangles. 
Since these objects are secondary to the main bubble-finding 
task, the site was designed so that they should be simple 
and quick to mark. Simple rectangles allow us to record the 
positions and approximate sizes of any interesting objects. 

Citizen scientists can discuss and share objects and im- 
ages via the 'Talk' interface 5 where more heavily discussed 
objects trend upwards as they do, for example, in a news 
aggregator. The 'Talk' web application is open source 6 and 
was developed by the Zooniverse team for general use on 
its projects, including the MWP. Through the use of 'Talk', 
interesting objects float to the top of discussion and are 
identified as interesting to the MWP scientists. Thus, by 
harnessing the social nature of the MWP, we can extract 
additional information from the classification process. 

Principally the objects highlighted by volunteers using 
'Talk' are visually interesting, unusual or defy classification 
in the primary interface. One example of such a feature are 
the 'yellowballs' - named by users because of their compact, 
rounded, yellow appearance in MWP images. It is believed 
that these represent a type of ultra-compact H n region and 
they are the subject of a future paper currently being pre- 
pared. They were an unintended consequence of our colour 
scheme and thus no flag was provided in the main interface 
to enable volunteers to mark their presence. 



2.1 Users' favourite bubbles 

Users are able to mark images as 'favourites' as they classify 
on the MWP. These images are often particularly beautiful, 
interesting or unusual. Ten of the most-favourited images are 
shown in Figure 3. The list includes many notable and well- 
known objects that have been 'rediscovered' by the MWP 
users, such as the Eagle Nebula (b) the Trifid Nebula (c), 
and the Galactic centre (h). Users not only select beautiful 



4 http://www.milkywayprojcct.org/tutorial 
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(g) Z=48.7°, Z>=-0.43°, zoom=3 (h) Z=359.58°, 6=-0.06°, zoom=2, towards Galactic centre 




(i) £=339.8°, 6=-0.13°, zoom=l, contains Westerlund 1 (j) i=51.8°, 6=0.63°, zoom=l 



Figure 3. Ten Milky Way Project images most-favourited by volunteers, in no particular order. Coordinates are image centres, image 
sizes are indicated by the zoom level (zoom). Zoom levels 1, 2 and 3 refer to images of 1.5 X 1°, 0.75 X 0.375° and 0.3 X 0.15° respectively. 
Colour figure available online. 
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Table 1. Tabic showing a user's score given the number of pre- 
cision bubbles they have drawn. Precision bubbles are those that 
require the use of multiple modifications to the default bubble 
parameters (see text for details). As such the number of precision 
bubbles drawn is a proxy for the care users are willing to take 
thus their experience with the tool. 
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images but also those that contain interesting objects. Im- 
age i contains the massive star cluster Westerlund 1 (West- 
erlund 1961), which is much-discussed on Milky Way Talk 
(see Section 2). 
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Figure 4. Flow chart describing the way that similar bubbles are 
found and combined. 



3 CATALOGUE CONSTRUCTION 

3.1 MWP user statistics and scoring 

By October 2011, over 35,000 people had logged in to milky- 
wayproject.org. 45% of those users classified at least one of 
the 12,263 MWP images; 25% classified 5 or more images 
and 5.7% classified 50 or more. Figure 2 shows the distribu- 
tion of users with the number of classifications they perform 
over the lifetime of their involvement with the project. 

The top-five nations visiting the Milky Way Project 
have been the United States (42%), United Kingdom (20%), 
Canada (4.3%), Poland (3.8%) and Germany (3.8%). The 
remaining 26.1% of project visitors have come from 173 
other countries. The MWP site is predominately accessed 
by English-speaking countries with a high level of internet 
connectivity, in addition it has been translated into Polish 
and hence that country also provides a large proportion of 
visitors. 

In order to assist in the data-reduction process, users 
are given scores according to how experienced they are at 
drawing bubbles. We treat the first 10 bubbles a user draws 
as practice drawings and these are not included in the fi- 
nal reduction. Users begin with a score of and are given 
scores according to the number of precision bubbles they 
have drawn (see Table 1). Precision bubbles are those drawn 
using the full toolset, meaning they have to have adjusted 
the ellipticity, the thickness and the rotation. This is done to 
ensure that users' scores reflect their ability to draw bubbles 
well. While only precision bubbles are used to score volun- 
teers, all bubbles drawn as included in the data reduction. 
The scores are used as weights when averaging the bubble 
drawings to produce the catalogue. 

3.2 Combining the user-drawn bubbles 

Combining all bubble drawings, at all zoom levels, the MWP 
has created a database of 520,120 user-drawn bubbles as 
of Oct 31st 2011. To identify bubbles with at least 5 user 
classifications for inclusion in our final catalogue, the data 



reduction process (shown as a flow chart in Figure 4) begins 
by splitting the dataset into 2° x 2° boxes and treating each 
box in turn. 

If a box contains 5 or more bubbles with a maximum 
outer ellipse that is between a half- and a whole-box, then 
a simple clustering algorithm picks out groups of these bub- 
bles with dispersions in their positions of less than a quar- 
ter of the box size (i.e. less than the radius of the smallest 
bubble under consideration). If a cluster contains at least 5 
bubble drawings, it is saved for additional processing and 
inclusion in the catalogue, and the bubble drawings are re- 
moved from the working list. Bubble drawings that are not 
clustered enough, or numerous enough, remain on the work- 
ing list for potential inclusion in a later iteration. The box 
is then split into four and the process repeats until no more 
boxes containing 5 or more bubbles are found, or until the 
box size falls below the smallest bubbles drawn on the MWP 
- the ellipse-drawing tool has a lower size limit of a diameter 
of 20 pixels (0.45' at the highest image zoom level). 

The same process is also run on an offset grid where 
the initial boxes are displaced by 1 degree in both galactic 
latitude and longitude. This catches bubbles that may fall 
on box boundaries. The two resultant lists of bubble-groups 
are combined later on by clustering bubbles that fall within 
0.5 radii separation from each other and which both have 
radii within 50% of each other. 

Each resulting cluster of bubbles marked > 5 times is 
combined into a single 'clean' bubble using a weighted mean, 
where the weighting is provided by the score of the user 
that drew each bubble (see Section 3.1). The bubble's mean 
size, position, angle (in degrees from North, in galactic co- 
ordinates) and thickness are all determined in this way (see 
Figure 5). The cleaned bubble catalogue is given in Table 4. 

The bubble's 'hit rate' is the ratio of the number of 
qualifying bubbles drawn to the number of times the bubble 
was seen by users on the MWP website. A clean bubble 
produced from a cluster of 5 user drawings placed onto an 
asset that was seen 50 times would have a hit rate of 0.1. 
In cases where a bubble could be marked in more than one 
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asset, for example across different zoom levels, the total view 
counts are summed, such that a cluster of 5 bubbles drawn 
onto two assets of 50 views each would have a hit rate of 
0.05. The hit rate gives a measure of consensus among users 
that a bubble is present in the data. 

The bubble's dispersion is also calculated as t he spread 
in coordinates of the individual classifications (•»/ a\ + erf, 
where a is the variance in the coordinate value). 

Only bubbles that were seen 50 times or more, and 
which have hit rates of 0.1 or more, are included in the final 
catalogue. This ensures that each final, cleaned bubble is a 
combination of at least 5 individual users' drawings and was 
drawn by volunteers at least 10% of the time when displayed 
on the website. 



3.3 Selection effects 

It is not known how many bubbles exist in the Galaxy, hence 
it is impossible to quantify the completeness of the MWP 
catalogue. There will be bubbles that are either not visi- 
ble in the data used on the MWP, or that are not seen as 
bubbles. Distant bubbles may be obscured by foreground 
extinction. Faint bubbles may be masked by bright Galac- 
tic background emission or confused with brighter nebular 
structures. Fragmented or highly distorted bubbles present 
at high inclination angles may not appear as bubbles to the 
observer. 

The MWP's 'citizen science' approach creates its own 
biases, whilst overcoming others experienced in similar stud- 
ies. By comparison with CP06, this study has many thou- 
sands of times more eyes scanning each section of the sky, 
and each section is broken down to an optimal colour stretch, 
thus improving the chances of seeing bubble- like structures. 
The majority of the MWP volunteers have no professional 
bias or expectation as to what constitutes a 'good' bubble. 
MWP volunteers may experience measurement fatigue when 
classifying assets with many bubbles. They may also suffer 
bright neighbour bias, and fail to draw quite obvious bubbles 
that are adjacent to very prominent or beautiful examples. 

3.4 Small bubbles and other objects 

In addition to marking elliptical bubbles on images, users are 
also encouraged to mark the locations of other interesting 
objects. Users can mark areas using a simple rectangle and 
are asked to label them as either a small bubble, green knot, 
dark nebula, star cluster, galaxy, fuzzy red object or other. 

The ellipse-drawing tool of the MWP has a lower size 
limit of a diameter of 20 pixels (0.45' at the highest image 
zoom level). The Small Bubble category allows users to mark 
bubbles which are too small to draw in detail but which can 
still be clearly made out. These small bubbles are reduced in 
a similar fashion to the more complex ellipses. To produce 
the catalogue of small bubbles listed in Table 3 we use only 
the small bubbles drawn by users at the highest zoom level 
(this is the vast majority of those drawn). Small bubbles 
marked at lower zoom levels are equivalent to larger bubbles 
at the higher zoom levels. By rounding their locations to the 
nearest 20 pixels, the drawings are clustered. Our catalogue 
of 1,362 small bubbles is given in Table 3. Each of these small 
bubbles was drawn by at least five users and was drawn by 




(a) Image: GLIMPSE only, as in CP06 



* 




(b) Image: GLIMPSE+MIPSGAL, from the MWP 




(c) Raw User Drawings ('Heat Map') 




(d) Reduced, 'cleaned' bubbles. 



Figure 5. Example of raw user drawings and reduced, cleaned 
result using a sample MWP image. A GLIMPSE-only colour sam- 
ple is included to illustrate the differences in the appearance of 
images inspected by CP06 and the MWP users. Image shown is 
centred at J=18.8°, 6=-0.125°, with size 1.5° X 0.75°. In image (c) 
all user-drawn bubbles are placed, from all zoom levels, with an 
opacity of 2.5%. In image (d) reduced bubbles are placed with an 
opacity 2x their hit rate, such that bubbles with hit rates > 50% 
are drawn as a solid white bubbles. Colour figure available online. 
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at least 10% of the volunteers who saw it - as with the main, 
large-bubble catalogue. 

Catalogues of green knots, dark nebulae, star clusters, 
galaxies, fuzzy red objects and objects in the other category 
are currently being prepared for later publication. 



4 RESULTS 

4.1 Catalogue description 

The final, reduced catalogue contains 5,106 visually identi- 
fied bubbles. These are split into a catalogue of 3,744 large 
bubbles drawn by users as ellipses, and a catalogue of 1,362 
small bubbles drawn by users at the highest zoom level im- 
ages in the MWP. These bubbles are plotted in Galactic 
coordinates in Figure 6. 

Each bubble in both lists has been drawn by at least five 
different individuals, and the listed parameters have been 
obtained from a weighted average based on each user's score 
(see Section 3.1). The complete catalogue can be accessed 
at http://data.milkywayproject.org. Table 4 gives the large 
bubbles, in order of hit rate. The large bubble catalogue 
includes values for the position in galactic longitude and 
latitude; mean geometric radius and mean thickness (as de- 
fined by CP06), position angle (given in degree from North), 
eccentricity (as defined by CP06); the hit rate (described in 
Section 3.2); and a hierarchy flag indicating whether bubbles 
are (1) identified as having further, smaller bubbles within 
their boundary, or (2) located on the rim of a larger bubble. 

4.2 Cross-Matching with Existing Catalogues 

For each bubble and small bubble produced by the MWP, 
cross-matching was performed with GLIMPSE bubbles (i.e. 
the CP06 and CWP07 catalogues) and the Paladini et al. 
(2003) and Anderson et al. (2011) catalogues of H II regions. 
Sources are marked as coincident when the central coordi- 
nate of the catalogue object lies within the radius of the 
MWP bubble. 12% of MWP bubbles matched GLIMPSE 
bubbles, and 86% of GLIMPSE bubbles were re-discovered 
by MWP. Similarly 10% and 7% of MWP bubbles are coin- 
cident with Paladini et al. (2003) and Anderson et al. (2011) 
H II regions, respectively. The MWP finds 86% of the avail- 
able Paladini sources and 96% of the Anderson sources. The 
presence of 24/im emission coincident with 20 cm emission 
was a selection criterion for Anderson et al. (2011) and so 
there may be some overlap in selection methods with the 
MWP. 

Mizuno et al. (2010) catalogued 416 disk and ringlike 
structures seen in the MIPSGAL 24 /im images and 
suggested that the majority of these objects were produced 
by evolved stars. The Mizuno et al. (2010) catalogue 
is dominated by small sources with radii <20", and so 
these objects should be unlikely to overlap with the 
main MWP bubbles catalogue. In fact we rediscover only 
9% of the sources in this catalogue and they constitute 
less than 1% of the combined MWP bubble catalogues. 
Of the 1,093 small MWP bubbles with a mean width 
of less than 1 arcmin, only five correspond to objects 
in the Mizuno et al. (2010) catalogue (for reference, 
the IDs of these bubbles are: MWP1G031730+07000S, 



Table 4. Crossover between the MWP and relevant catalogues of 
bubbles (CP06; CWP07), H II regions (Paladini et al. 2003; An- 
derson ct al. 2011), and MIPSGAL ring-like structures (Mizuno 
et al. 2010). 



Catalogue 


Fraction MWP 
Bubbles Matched 


Fraction 
Rediscovered 


CP06 and CWP07 


0.12 


0.85 


Paladini ct al. (2003) 


0.10 


0.86 


Anderson ct al. (2011) 


0.07 


0.96 


Mizuno ct al. (2010) 


0.01 


0.09 



MWP1G314360+04900S, MWP1G319220+01600S, 
MWP1G334110+03800S, MWP1G358770+01100S). 

The above crossover fractions (summarised in Table 4) 
show that the MWP has excellent overlap with existing bub- 
ble catalogues and is also more complete, in terms of locating 
H II regions, than the two Churchwell studies. The lack of 
agreement between the MWP and Mizuno et al. (2010) is ev- 
idence that the small bubble catalogue is not contaminated 
by small, 24 fj,m disk- and ring-like structures. 

4.3 Errors 

A number of quality control measures, such as the user 
weighting scheme, were adopted from the outset; these 
are described in Section 3.1. To assess the performance of 
the processing procedure, we examine the four bubbles in 
the catalog with the highest classification scores, i.e. those 
drawn by the largest number of users. These bubbles were 
identified by an average of 243.2 users. Three of these were 
previously identified as bubbles and HII regions by other au- 
thors (CP06; Misanovic et al (2002); Lockman et al (1996); 
Kuchar & Clark (1997)). Bubble MWP1G303056+01645 
(Figure 10) has no associated CP06 bubble. 

In bubbles MWP1G309059+01661 (Figure 7) and 
MWP1G303056+01645 (Figure 10) the computed weighted 
averages for the bubble position and size parameters, in- 
dicated with the dashed lines, accurately follow the distri- 
bution of individual classifications, and the distribution is 
well approximated by a gaussian distribution. In these cases, 
each user classification can be reasonably considered to be 
an independent measurement of the same quantity. In other 
cases, such as with MWP1G304463-00217 (Figure 9), there 
are skews or multiple peaks in the distribution of classifica- 
tions, indicating a relative lack of consensus among the users 
over the location, size, shape or multiplicity of the bubble. 
Indeed, multiple peaks in the distribution may suggest the 
presence of more than one bubble at a given location. 

These variations can be attributed to the definition of 
our clustering algorithm. Using a tighter clustering thresh- 
old will do better at separating out closely spaced bubbles, 
but will also artificially fragment single bubbles into multi- 
ples where the bubble's positional coordinates have a high 
uncertainty value. Our data reduction could be modified to 
track the level of dispersion in bubble-drawing clusters and 
dynamically split or reject clusters based on this value. We 
aim to address this issue in a subsequent data release. 

To assess the reliability of the data presented in a quan- 
titative way, two reliability metrics were defined: 
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Table 2. Table of bubble parameters (shown here are the 50 bubbles from the catalogue with the highest hit rates). The complete 
large bubble catalogue contains 3,744 visually identified bubbles. Where cross-correlation is possible, identifiers from CP06 and CWP07 
are given. Hierarchy flags denote (1) bubbles identified as having smaller bubbles on their rim, and (2) bubbles located within a larger 
bubble. The complete bubble catalogue can be found online at http://data.milkywayproject.org. 
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Figure 6. All 5,106 MWP bubbles plotted in Galactic coordinnates. The MWP large bubble catalogue is marked with outer radii as 
grey ellipses, small bubbles shown as green crosses. 
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Table 3. Tabic of MWP small bubbles. Shown here are the 50 bubbles from the small bubble catalogue with the highest hit rates. The 
complete small bubble catalogue contains 1,362 visually identified bubbles. Where cross-correlation is possible, identifiers from CP06 and 
CWP07 are given. Hierarchy flags denote (1) bubbles identified as having smaller bubbles on their rim, and (2) bubbles located on the 
rim of a larger bubble. The complete bubble catalogue can be found online at http://data.milkywayprojcct.org. 
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(i) Hit rate: the fraction of all users that were presented 
with the region of sky, who drew a bubble at this location 
(described in Section 3.2). This is a measure of the reliability 
of the existence of a bubble; and 

(ii) Dispersion: the spread in coordinates of the individ- 
ual classifications ( \Jo\ + erf , where a is the variance in the 



coordinate value). This describes the uncertainty on a bub- 
ble's location. 

Note that these metrics do not offer any insight into the 
physical nature of the bubbles - i.e. how likely a feature with 
the form of a bubble is to actually be an H II region - rather 
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Figure 7. Error measurements for MWP bubble MWP1G309059+01661. This bubble has a hit rate of 0.437, and a dispersion of 1.61'. 
Top figures show reduced and raw bubble drawings. Bottom figures show dispersions in measurements of position and size. Colour figure 
available online. 



they reflect the level of consensus gathered from the users 
that a bubble appears at this location or with this form. 

The dispersion measure shows the varying levels of un- 
certainty on the bubble positions gathered from the classifi- 
cations. In addition, the clustering threshold chosen in our 
data processing algorithm may merge bubbles that closer 
inspection suggest to consist of several bubble components. 
Thus the dispersion metric offers one way of assessing how 
many bubbles may have been 'lost' in data processing. 

Figure 11 shows the cumulative distribution function 
of the dispersion of all large bubbles, normalised to their 
effective radii. If we posit that a dispersion on a bubble's 
central position covering much of the bubble's effective ra- 
dius is highly suggestive of multiplicity, we can estimate how 
many bubbles may be lost. While 89% of large bubbles have 
a dispersion < 0.5 R e //, 2.4% (88 bubbles) show a disper- 
sion > 0.75 R e //- Thus assuming that each of these is highly 
likely to host at least one additional bubble around its rim, 
we can estimate that around 80-100 bubbles were artificially 
merged by the algorithm. 

4.4 Bubble catalogue properties 

The longitudinal distribution of the MWP bubbles is shown 
in Figure 12. This figure shows a broad rise and fall either 
side of the Galactic centre, with the number of bubbles be- 
ginning low near the Galactic centre, rising, and then dimin- 
ishing toward the edge of the survey at I ± 65° . The notable 
lack of bubbles around I = 60° is partly due to reaching the 



survey's edge at I — 65° and also due to an absence of any 
filaments or bubbles around I = 58°. Figure 12 marks sev- 
eral notable Galactic line-of-sight features in red, adopted 
from Beuther et al. (2011). Many of the rises and falls in 
the number of bubbles across the range of Galactic longi- 
tude appear to derive from the large-scale structure of the 
Milky Way. 

The drop in bubble count near the Galactic centre may 
be a physical effect but could also arise due to confusion from 
background emission toward the centre of the Milky Way. 
Beuther et al. (2011) studied the distribution of submm 
clumps from ATLASGAL - The APEX Telescope Large 
Area Survey of the GALaxy (Schuller et al. 2009). They find 
a large peak in the number of sources toward the Galactic 
centre and note that this is in contrast with current surveys 
of H II regions (e.g. Anderson et al. 2011) and recent surveys 
of H 2 and CH3OH masers (Green et al. 2011; Walsh et al. 
2011). 

The distribution of large bubbles with latitude is shown 
in Figure 13. We divide the GLIMPSE/MIPSGAL survey 
area by longitude into northern (I = 0° to 65°) and south- 
ern (I = 295° to 360°) regions to facilitate comparison of 
possible morphological, positional, and size differences with 
CP06 and CWP07. The same profile is shown in CP06 Fig- 
ure 5 and the plots share the same asymmetry toward lower 
latitudes. 

Figures 14, 15 and 16 show the distributions of the 
MWP bubble radii, thicknesses and eccentricities. Compar- 
ing these properties with those derived in CP06 reveals some 
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Figure 8. Errors for bubble MWP1G050955+01074. This bubble has a hit rate of 0.282, and a dispersion of 4.31'. See Figure 7 for more 
information. Colour figure available online. 
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Figure 9. Errors for bubble MWP1G304463-00217. This bubble has a hit rate of 0.395, and a dispersion of 0.70'. See Figure 7 for more 
information. Colour figure available online 
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Figure 10. Errors for bubble MWP1G303056+01645. This bubble has a hit rate of 0.355, and a dispersion of 0.59'. See Figure 7 for 
more information. Colour figure available online. 
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Figure 11. Cumulative distribution function of all large MWP 
bubbles normalised to their effective radii. 



differences. The MWP bubbles show a greater range in ra- 
dius and thickness, though both plots show broadly simi- 
lar features to the equivalents in CP06. Most strikingly, the 
MWP eccentricities peak at ~ 0.35, compared to a value of 
~ 0.65 in CP06. This could be a result of averaging the 
parameters of multiple ellipses to create the MWP large 
bubbles, i.e. merging a large number of elliptic annuli tends 
toward circularity as the number and variety increases. 

Figure 17 shows the distribution of the ratio of bubble 
thickness to outer diameter. Compared to CP06 Figure 10, 
the MWP catalogue displays a wider range of values, and 



in general the MWP bubbles are slightly thicker relative to 
their size (peak value at ~ 0.4) than those in CP06 (peak 
value at ~ 0.25). This may be partially due to the addition 
of 24 fim MIPSGAL data in the MWP, which often gives 
bubbles a 'fuzzier' rim - see Figure 5a and b. 



4.5 Bubble Distances 

Cross-matching with Anderson & Bania (2009) provides dis- 
tances to 185 of the larger MWP bubbles (cross-matching as 
described in Section 4.2). Figure 18 shows these 185 bubbles 
plotted in terms of their diameters against (a) their distance 
from the Sun and (b) their distance from the Galactic cen- 
tre. Figure 18a shows a slight tendency for more distant 
bubbles to be larger - most likely a reflection of the selec- 
tion effect whereby only very large distant bubbles are easily 
seen in the GLIMPSE/MIPSGAL images used in the MWP. 
Figure 18b shows little correlation other than reflecting the 
fact that fewer bubbles are seen at greater distances from 
the Galactic centre - this is to be expected given the longi- 
tudinal range of the MWP and the confusion effect where 
nearer bubbles and dust obscure more distant ones, looking 
toward the Galactic centre. 

Figure 19 shows that larger bubbles tend to have thin- 
ner shells relative to their diameters. This could be the effect 
of material cooling and condensing as the bubble expands. 
The apparent relationship warrants further investigation but 
we note caution that this relationship may be biased by the 
effect of a minimum size and thickness of a bubble in the 
MWP drawing tool. 
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Figure 12. Histogram showing distribution of combined MWP small- and large-bubble catalogues with Galactic longitude I. Notable 
line-of-site features (including Galactic spiral arms) marked as dashed boxes. 



4.6 'Heat maps' 

In addition to the reduced bubble catalogue, a crowd- 
sourced 'heat map' of bubble drawings has also been pro- 
duced. This simple map reflects the full range of classi- 
fications placed onto the MWP images. All 520,120 bub- 
bles drawn by all users are placed onto the sky with 
an opacity of 2.5% meaning that 40 individuals need to 
have drawn over the same region for it to become fully 
opaque (white). Examples of these images are shown in Fig- 
ures 5 and 20. These MWP 'heat maps' are available at 
http://data.milkywayproject.org as FITS and DS9 region 
files. 

The MWP 'heat maps' allow the bubble drawings to be 
explored without them needing to be reduced to elliptical 
annuli. Rather, the 'heat maps' allow contours of overlap- 
ping classifications to be drawn over regions of the Galac- 
tic plane reflecting levels of agreement between independent 
classifiers. In most cases the structures outlined in these 
maps are photo-dissociation regions traced by 8 (im emis- 
sion, but more fundamentally they are regions that multiple 
volunteers agree reflect the rims of bubbles. 




] -0.5 0.5 
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Figure 13. Distribution of MWP large bubbles with b (see figure 
5 of CP06). In additional to the overall distribution (drawn as a 
solid line) the dashed line shows just southern bubbles (I = 295° 
to 360°) and the dotted line shown just northern bubbles (Z = 0° 
to 65°). 



5 DISCUSSION 

CP06 and CWP07 performed visual classification, using a 
handful of experts, of GLIMPSE data only. The MWP im- 
ages include data from the MIPSGAL 24 /im survey, which 
can enhance the shape and definition of bubble-like struc- 
tures (see Figure 5). Most likely the biggest difference be- 
tween the MWP and GLIMPSE studies is that tens of thou- 
sands of classifiers are involved in the MWP and most begin 
with no grounding in the astrophysical phenomena they are 
told to locate. This results in a set of unbiased classifications, 
sourced from multiple independent classifiers who are able 
to explore and annotate the entire dataset without fatigue. 



5.1 The nature of MWP bubbles 

Bubbles were identified in the MWP based on ringlike or 
arclike shapes observed by multiple independent classifiers 
in infrared images of the Galactic plane. Many astrophysical 
phenomena can give rise to such features. As CP06 noted, 
aside from young, massive stars several classes of evolved ob- 
jects can produce bubbles in the interstellar medium (ISM), 
including supernova remnants (SNRs), planetary nebulae 
(PNe), and asymptotic giant branch (AGB) stars. SNRs and 
AGB stars don't excite PAH emission, thus these objects 
generally lack the 8 /im shells characteristic of the GLIMPSE 
bubbles. Regardless, 3 SNRs were included in the CP06 cat- 
alogue and a few bright, round 24 (im SNRs may be included 
in the MWP catalogue. 

Other lines of evidence also suggest that the MWP bub- 
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Figure 14. Distribution of angular diameter of MWP large bub- 
bles - values from figure 6 of CP06, scaled to 10 X, shown as 
dashed line. The smallest radius users can draw is 0.45'. The 
largest image served to users is 1.5 X 1°. 
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Figure 15. Distribution of thicknesses of MWP large bubbles - 
values from figure 9 of CP06, scaled to lOx , shown as dashed line 
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Figure 16. Distribution of eccentricities of MWP large bubbles 
- values from figure 13 of CP06, scaled to lOx, shown as dashed 
line. 



Figure 17. Distribution of ratio of thickness to outer diameter 
for MWP large bubbles - values from figure 10 of CP06, scaled 
to lOx, shown as dashed line. 

bles, like the GLIMPSE bubbles, predominantly trace mas- 
sive star formation. The distribution of bubbles with Galac- 
tic latitude reflects a low scale height (Figure 13), similar 
to molecular clouds and the known Galactic OB population 
(see CP06). The CP06 bubble catalogue contained only 12% 
of available Paladini et al. (2003) H II regions. The MWP 
catalogue contains 86% of the Paladini sources, indicating 
that the new catalogue is more complete. 

A bubble is produced around a massive star when an 
H II region, driven by thermal overpressure, stellar winds, ra- 
diation pressure, or a combination of these feedback mech- 
anisms, expands into the surrounding cold ISM, sweeping 
up gas and dust into a dense shell surrounding a low- 
density, evacuated cavity (Weaver et al. 1977; Garcia-Segura 
& Franco 1996; Arthur et al. 2011; Draine 2011). The rela- 
tive contributions of different feedback mechanisms likely de- 
pend on the properties of the driving star(s), with the most 
massive, early O-type stars combining powerful stellar winds 
with high UV luminosities producing 'wind-blown bubbles,' 
while lower-mass, late-O and B dwarfs give rise to 'classi- 
cal' H II regions powered by UV photons alone (CWP07; 
Watson et al. 2008). Castor et al. (1975) and Weaver et al. 
(1977) derived analytic solutions for the expansion of stel- 
lar wind-blown bubbles into a uniform low-density medium. 
More recent modelling efforts have included the effects of 
ionising radiation and the stellar winds (e.g. Capriotti & 
Kozminski 2001; Draine 2011). 

The wind-blown bubbles around massive stars produced 
by these models display the following general structure: 

• An inner cavity cleared rapidly by a freely flowing hy- 
personic stellar wind; 

• A high-temperature region of shocked stellar wind ma- 
terial (T > 10 6 K); 

• A shell of shocked, photo-ionised gas (T ~ 10 4 K); 

• An shell of non-shocked, ionised gas (T ~ 10 4 K). 

• An outer shell of neutral material. 

The bright PAH emission in the PDRs surrounding 
H II regions produces the bright 8 /tm bubble rims, while 
dust mixed with the ionised gas and heated by the hard ra- 
diation field produces 24 jj,m emission interior to the bubbles 
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Figure 18. Bubble diameters against (a) distance from the Sun and (b) distance from the Galactic centre, for 185 MWP bubbles 
cross-matched with Anderson & Bania (2009). 



(Povich et al. 2007; Watson et al. 2008; Everett & Church- 
well 2010; Watson et al. 2010). These considerations guided 
our choice of multiband colour combination for the MWP 
images, which visually bias identification toward H n re- 
gions. This is consistent with the cross-matching results in 
Section 4.2. 

Watson et al. (2008) proposed that wind-blown bubbles 
with cleared, central cavities tend to produce a toroidal or 
arc-shaped morphology in the 24 /im emission, while this 
emission is more likely to be centrally peaked in classical 
H II regions. Draine (2011) found that radiation pressure 
alone can produce cleared central cavities in H II regions, 
but suggested cases where winds must also play a role. A 
wide range of 24 jj,m morphologies is apparent among the 
MWP bubbles. 

The models described above consider the effects of a sin- 
gle star on its surrounding uniform-density medium. Mas- 
sive stars form preferentially in clusters within molecular 



clouds with highly non-uniform densities. Freyer et al. (2003) 
performed two-dimensional radiation hydrodynamics simu- 
lations to study the combined influence of a massive star's 
stellar wind and its ionising radiation on the evolution of 
the circumstcllar cloud material. Mellema et al. (2006) sim- 
ulated the evolution of an H II region surrounding a single 
O star, focusing on the interaction of the ionisation front 
with the turbulent cloud medium. More recently, Arthur 
et al. (2011) performed magnetohydrodynamic (MHD) sim- 
ulations of H II regions around OB stars expanding into 
a turbulent, magnetised medium. All of these simulations 
show strong inhomogeneities around the boundary of the 
expanding shell, reproducing the clumps, globules and fila- 
ments observed in H II regions at optical or IR wavelengths 
and seen in the MWP bubbles. 
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Figure 19. Variation of the thickness-diameter ratio of large 
MWP bubbles (Table 4) with their physical size, for 185 MWP 
bubbles cross-matched with Anderson & Bania (2009). 



5.2 Multiple bubbles and potentially triggered 
star formation 

A major application of the CP06 and CPW07 bubble cata- 
logs has been triggering studies (e.g. Deharveng et al. 2008; 
Watson et al. 2008, 2010; Zavagno et al. 2010, and refer- 
ences therein), as the Churchwell catalogs provided a large 
sample of (mostly previously unknown) regions of potential 
triggering by massive stars. In any individual source, diffi- 
culties with establishing cause-and-effect are significant, and 
the Churchwell catalogs provided a basis for triggering inves- 
tigations to move beyond case studies of individual objects 
to statistics. 

Bubbles serve as laboratories to test theories of sequen- 
tial, massive star formation triggered by massive star winds 
and radiation pressure (Elmegreen & Lada 1977; Whitworth 
et al. 1994). CP06 and CWP07 noted that a (small) fraction 
of bubbles exhibited hierarchical structure, meaning that 
one or more small, 'daughter' bubbles were found on the 
rims of, or projected inside, larger, 'parent' bubbles. The 
prevalence of triggering is a key unresolved question in the 
study of massive star formation, with important implica- 
tions for extragalactic studies as well as detailed star forma- 
tion physics. Quantifying triggering allows the formulation 
of galaxy-scale star formation rate prescriptions for use in 
simulations, particularly if the physical conditions in trig- 
gered regions affect the IMF, as suggested by Whitworth 
et al. (1994) and Dale et al. (2009). 

Rahman & Murray (2010), however, observed that 
Galactic giant H II regions were arranged on the rims of 
very large (up to 100 pc diameter), mid-IR bubble structures 
associated with the most luminous sources of free-free emis- 
sion in the Galaxy, suggesting large-scale triggering driven 
by very massive clusters and OB associations. 

However, whether these structural features are a direct 
result of the expanding shell's interaction with the molecular 
cloud, or whether the clearing of material simply reveals the 
underlying turbulent cloud structure, is not yet clear. Beau- 
mont & Williams (2010) proposed that the 3-dimensional 
geometry of bubbles in the CP06 catalogue, when traced by 
molecular gas, resembles flattened rings rather than spher- 



ical shells, and they noted that this geometry could reduce 
the efficiency of triggering. 

From results of smoothed particle hydrodynamic (SPH) 
simulations of the ionising radiation from massive stars or 
small protoclusters, Dale & Bonnell (2011) argue that ISM 
bubbles are features of the turbulent nature of the molecular 
clouds rather than shells created by feedback from massive 
stars and clusters. Their simulation involves multiple clus- 
ters within a very massive molecular cloud. In this regard 
their simulation is more analogous to a region such as the 
Carina Nebula than to the kind of bubbles predominantly 
seen in the MWP, where the ionisation is dominated by one 
or two stars. 

Most studies to date have asked the question, 'what 
percentage of bubbles show evidence for triggered star for- 
mation?' Understanding Galactic-scale star formation in the 
Milky Way, however, requires answering not this question 
but its corollary: 'what fraction of all (massive) star for- 
mation is triggered?'. Thompson et al. (2011) have recently 
made the first effort to address this question, using the Red 
MSX Source (RMS) database of massive YSOs and the bub- 
bles from CP06 and CWP07, but the major caveat of their 
analysis is that the Churchwell bubble catalogs are incom- 
plete. By identifying an order of magnitude more bubbles, 
and providing a reliability indicator in the form of the hit 
rate, the MWP bubble catalog greatly ameliorates this prob- 
lem. 

In addition, Oey et al. (2005) have suggested that the 
most convincing candidates for triggered star formation are 
regions where three-generation hierarchies can be estab- 
lished. The MWP has identified a much larger number of 
broken, old bubbles and small bubbles than the CP06 and 
CWP07 catalogs. As such, it is an excellent dataset for 
searching for multiple hierarchies. Indeed, 29% of MWP 
bubbles have hierarchy flags of 1 or 2, indicating bubbles 
that are located on or within larger bubbles and bubbles 
that have smaller bubbles situated on or within them. 

Additional analysis to assess triggering related to MWP 
bubbles (using multiwavelength datasets) is ongoing and will 
be the subject of a subsequent paper. 



6 SUMMARY AND FUTURE WORK 

A new catalogue of 5,106 infrared bubbles has been created 
through visual classification via the MWP website. Bubbles 
in the new catalogue have been independently measured by 
at least 5 individuals, producing averaged parameters for 
their position, radius, thickness, eccentricity and position 
angle. Citizen scientists have independently rediscovered the 
locations of 86% of the Churchwell et al. (2006) and Church- 
well et al. (2007) bubble catalogues and 96% of the Anderson 
ct al. (2011) H II region catalogue, whilst finding an order 
of magnitude more objects. 

The MWP bubble catalog constitutes a resource that, in 
combination with other recent and ongoing Galactic Plane 
surveys of star-formation tracers (including the Methanol 
Multibeam Survey (MMB), the Bolocam Galactic Plane 
Survey (BGPS), ATLASGAL, and HiGAL), has the poten- 
tial to provide sufficient statistics to address the question 
of how prevalent, and important, triggered star formation 
really is. In addition we hope that this new resource will 
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(e) Z=31.8°, ft=-0.06°, zoom=2 (Heat Map) 
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Figure 20. Examples of sites showing potential evidence of triggering. In each case the 'heat map' and reduced data is shown overlaid 
on the MWP image. Coordinates are image centres, image sizes are indicated by the zoom level (zoom). Zoom levels 1, 2 and 3 refer to 
images of 1.5 X 1°, 0.75 X 0.375° and 0.3 X 0.15° respectively. Colour figure available online. 



complement these surveys as a tracer of massive star forma- 
tion on Galactic scales. 

Also outlined is the creation of a 'heat map' of star- 
formation activity in the Galactic plane. This online resource 
provides a crowd-sourced map of bubbles and arcs in the 
Milky Way, and should enable better statistical analysis of 
nearby star- format ion sites. 

Additonal papers are currently being prepared to out- 
line catalogues of 'green knots', dark nebulae, star clusters, 
galaxies and 'fuzzy red objects' that have also been created 
by the MWP's community of citizens scientists. Similarly, we 
anticipate a second, refined bubble catalogue incorporating 
not only better data-reduction techniques but also 100,000s 
more bubble drawings by volunteers. 



ACKNOWLEDGMENTS 

This publication has been made possible by the partici- 
pation of more than 35,000 volunteers on the Milky Way 
Project. Their contributions are acknowledged individually 



at http://www.milkywayproject.org/authors. We would like 
to thank Bob Benjamin for his helpful suggestions. 

The Milky Way Project, and R.J.S. were supported by 
The Lever hulme Trust. The 'Talk' discussion tool used in 
the MWP was developed at the Adler Planetarium with 
support from the National Science Foundation CDI grant : 
DRL-0941610. M.S. P. was supported by a National Science 
Foundation Astronomy & Astrophysics Postdoctoral Fellow- 
ship under award AST-0901646. 

C.J.C. is supported by an NSF Astronomy and Astro- 
physics Postdoctoral Fellowship under award AST-1003134. 
Support for the work of K.S. was provided by NASA through 
Einstein Postdoctoral Fellowship grant number PF9-00069 
issued by the Chandra X-ray Observatory Center, which is 
operated by the Smithsonian Astrophysical Observatory for 
and on behalf of NASA under contract NAS8-03060. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. 



© 0000 RAS, MNRAS 000, 000-000 



Milky Way Project DR1 19 



REFERENCES 

Anderson, L. D., & Bania, T. M. 2009, ApJ, 690, 706 
Anderson, L. D., Bania, T. M., Balser, D. S., & Rood, R. T. 

2011, VizieR Online Data Catalog, 2194, 40032 
Arthur, S. J., Henney, W. J., Mellema, G., de Colle, F., & 

Vazquez-Semadeni, E. 2011, MNRAS, 414, 1747 
Bania, T. M., Anderson, L. D., Balser, D. S., & Rood, R. T. 

2010, ApJL, 718, L106 

Beaumont, C. N., & Williams, J. P. 2010, ApJ, 709, 791 
Beaumont, C. N., Williams, J. P., & Goodman, A. A. 2011, 

ArXiv e-prints 
Bendo, G. J., et al. 2008, MNRAS, 389, 629 
Benjamin, R. A., et al. 2003, PASP, 115, 953 
Beuther, H., et al. 2011, ArXiv e-prints 
Capriotti, E. R., & Kozminski, J. F. 2001, PASP, 113, 677 
Cardamone, C, et al. 2009, MNRAS, 399, 1191 
Carey, S. J., et al. 2009, PASP, 121, 76 
Castor, J., McCray, R., & Weaver, R. 1975, ApJL, 200, 

L107 

Churchwell, E., et al. 2006, ApJ, 649, 759 (CP06) 
— . 2007, ApJ, 670, 428 (CWP07) 
— . 2009, PASP, 121, 213 

Dale, J. E., & Bonnell, I. 2011, MNRAS, 414, 321 

Dale, J. E., Wiinsch, R., Whitworth, A., & Palous, J. 2009, 

MNRAS, 398, 1537 
Deharveng, L., Lefloch, B., Kurtz, S., Nadeau, D., Pomares, 

M., Caplan, J., & Zavagno, A. 2008, A&A, 482, 585 
Draine, B. T. 2011, ApJ, 732, 100 
Elmegreen, B. G., & Lada, C. J. 1977, ApJ, 214, 725 
Everett, J. E., & Churchwell, E. 2010, ApJ, 713, 592 
Flagey, N., Boulanger, F., Noriega-Crespo, A., Paladini, 

R., Montmerle, T., Carey, S. J., Gagne, M., & Shenoy, S. 

2011, A&A, 531, A51+ 

Freyer, T., Hensler, G., & Yorke, H. W. 2003, ApJ, 594, 
888 

Garcia-Segura, G., & Franco, J. 1996, ApJ, 469, 171 

Green, J. A., et al. 2011, ApJ, 733, 27 

Joy, K., et al. 2011, Astronomy and Geophysics, 52, 020000 

Lintott, C, et al. 2011, MNRAS, 410, 166 

Lintott, C. J., et al. 2008, MNRAS, 389, 1179 

— . 2009, MNRAS, 399, 129 

Mellema, G., Arthur, S. J., Henney, W. J., Iliev, I. T., & 

Shapiro, P. R. 2006, ApJ, 647, 397 
Mizuno, D. R., et al. 2010, AJ, 139, 1542 
Oey, M. S., Watson, A. M., Kern, K., & Walth, G. L. 2005, 

AJ, 129, 393 

Paladini, R., Burigana, C, Davies, R. D., Maino, D., 
Bersanelli, M., Cappellini, B., Platania, P., &: Smoot, G. 
2003, A&A, 397, 213 

Povich, M. S., et al. 2007, ApJ, 660, 346 

— . 2009, ApJ, 696, 1278 

Raddick, M. J., Bracey, G., Gay, P. L., Lintott, C. J., Mur- 
ray, P., Schawinski, K., Szalay, A. S., & Vandenberg, J. 
2010, Astronomy Education Review, 9, 010103 
Rahman, M., & Murray, N. 2010, ApJ, 719, 1104 
Relano, M., & Kennicutt, Jr., R. C. 2009, ApJ, 699, 1125 
Schuller, F., et al. 2009, A&A, 504, 415 
Smith, A. M., et al. 2011, MNRAS, 412, 1309 
Thompson, M. A., Urquhart, J. S., Moore, T. J. T., & 

Morgan, L. K. 2011, ArXiv e-prints 
Walsh, A. J., et al. 2011, MNRAS, 416, 1764 



Watson, C, Hanspal, U., & Mengistu, A. 2010, ApJ, 716, 
1478 

Watson, C, et al. 2008, ApJ, 681, 1341 

Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, 

R. 1977, ApJ, 218, 377 
Westerlund, B. 1961, AJ, 66, 57 

Whitworth, A. P., Bhattal, A. S., Chapman, S. J., Disney, 
M. J., & Turner, J. A. 1994, MNRAS, 268, 291 
Zavagno, A., et al. 2010, A&A, 518, L101 



© 0000 RAS, MNRAS 000, 000-000 



